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ABSTRACT 

A  cantilever  beam  force  transducer  was  modelled  as  a  massless 
elastic  section  at  the  base  with  the  remaining  section  of  the  beam 
rigid  and  having  mass.  Computer  programs  were  written  to  simulate 
free  or  forced,  damped  or  undamped  vibrations  of  the  beam.  Good  agree¬ 
ment  was  found  between  predicted  and  experimental  frequencies  of  un¬ 
damped  free  vibration  for  two  different  beams.  After  further  veri¬ 
fication,  the  computer  programs  can  be  used  to  determine  beam  configu¬ 
rations,  viscous  damping  factors,  and  loading  rates  which  will  reduce 
unwanted  oscillations  of  the  transducer  element. 
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1.  INTRODUCTION 

The  Chemistry  Saction  at  the  Defence  Research  Establishment 
Suffield  is  studying  the  viscoelastic  properties  of  fluids.  It  was 
noticed  that  a  pencil  point  quickly  lifted  from  a  fluid  surface  would 
pull  up  a  filament  the  size  of  which  was  an  indication  of  the  fluid's 
viscoelastic  properties.  Based  on  that  observation,  two  members  of  the 
technical  staff  at  DRES  built  the  measurement  system  shown  schematically 
in  Figure  1 . 

The  pencil  point  was  replaced  by  a  short  length  of  .055-inch 
diameter  steel  rod  called  the  T-piece.  Instead  of  lifting  the  T-piece 
from  the  fluid  surface,  the  fluid  contained  in  a  1-inch  diameter  by 
0.1-inch  deep  cup  was  lowered  away  from  the  T-piece.  The  cup  rested 
upon  the  plunger  of  a  syringe.  Rate  of  descent  of  the  plunger  was  con¬ 
trolled  by  the  metering  valve  and  vacuum  pump. 

The  weight  of  the  fluid  filament  applied  a  vertical  load  to  the 
tip  of  the  cantilever  beam  causing  strain  in  the  beam  which  in  turn  caused 
a  change  in  electrical  resistance  of  the  strain  gauges  proportional  to 
filament  weight.  As  the  filament  would  reach  a  maximum  weight  and  then 
begin  to  drip  back  into  the  cup,  so  the  strain  versus  time  signal  would 
show  a  sharp  rise  followed  by  a  gradual  decline.  However  the  loading 
function  excited  oscillations  of  the  beam  which  tended  to  obscure  the 
desired  signal. 
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Several  remedies  were  considered.  A  small  amount  of  grease 
placed  between  the  beam  and  one  of  its  support  members  was  an  effective 
damper,  but  the  grease  slowly  oozed  out  and  it  was  difficult  to  maintain  the 
same  degree  of  damping  from  day  to  day.  First  attempts  at  electronically 
filtering  the  signal  had  results  similar  to  mechanical  overdamping:  the 
unwanted  oscillations  were  removed,  but  so  was  most  of  the  signal.  A  Mark  1 
beam  and  a  Mark  2  beam  of  slightly  different  dimensions  were  tried.  While 
both  had  an  oscillation  problem,  the  Mark  2  beam  vibrated  at  a  higher  fre¬ 
quency  and  with  greater  amplitude.  Hence  an  analysis  was  attempted  in  order 
to  determine  what  beam  dimensions  would  minimize  the  oscillation  problem. 

The  author's  approach  was  to  predict  the  behaviour  of  beam  models 
with  the  IBM  Continuous  System  Modeling  Program  (CSMP).  Figure  2  shows  the 
Mark  1  beam  and  its  model.  Elastic  deformation  was  assumed  to  occur  in  the 
section  of  beam  closest  to  the  support.  The  mass  of  this  section  was  neg¬ 
lected  because  deflections,  accelerations  and  hence  inertia  forces  were  low 
compared  to  those  of  the  other  part  of  the  beam.  The  flanged  section  of  beam 
was  assumed  rigid  and  its  mass  significant  because  of  relatively  large  de¬ 
flections  and  inertia  forces.  In  calculating  the  mass  of  the  rigid  section, 
the  masses  of  any  attached  damper,  the  T-piece  and  fluid  filament  were  neg¬ 
lected.  Some  compensation  for  those  assumptions  was  made  by  ignoring  the 
absence  of  mass  in  lightening  holes.  Table  I  defines  the  beam  nomenclature, 
and  Figure  3  illustrates  the  coordinate  system.  Because  the  position  of  the 
rigid  section  could  be  specified  by  two  coordinates,  the  model  was  a  two- 
degree-of-freedom  system  and  two  natural  frequencies  of  vibration  were  expected. 

This  report  describes  development  of  a  CSMP  "Vibrating  Cantilever" 
program  and  an  auxiliary  FORTRAN  program.  Instructions  for  their  use  are 
given.  Users  unfamiliar  with  CSMP  will  find  additional  background  information 
in  Reference  1.  Free  or  forced  vibration,  with  or  without  damping,  can  be 
examined.  Good  agreement  was  found  between  predicted  and  experimental  fre¬ 
quencies  of  undamped  free  vibration  for  the  Mark  1  and  Mark  2  beams.  Further 
verification  of  the  program  is  suggested  in  Section  7.  The  most  desirable 
beam  configuration,  viscous  damping  factor,  and  loading  rate  can  then  be 
found  by  examining  the  properties  of  various  candidate  configurations. 
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TABLE  I 


BEAM  NOMENCLATURE 


Hi  =  length  of  elastic  portion  of  beam,  in 

i2  =  length  of  clamp  to  centre  of  gravity  of  rigid  portion,  in 

H 3  =  length  of  clamp  to  line  of  action  of  damping  force,  in 

st4  =  length  of  clamp  to  line  of  action  of  forcing  function,  in 

b!  _  width  of  elastic  portion,  in 

b2  =  width  of  rigid  portion,  in 

dj  =  thickness  of  elastic  portion,  in 

d2  =  thickness  of  rigid  portion,  in 

W  =  specific  weight  of  rigid  material,  lb^  in-3 

E  =  Young's  modulus  of  elastic  material,  lb^  in-2 

2.  DERIVATION  OF  EQUATIONS  OF  MOTION 

Free  body  diagrams  of  the  elastic  and  rigid  portions  of  the 
beam  are  given  in  Figure  4.  It  can  be  shown  that  an  elastic  cantilever 
subjected  to  a  vertical  load  P  and  bending  moment  M  at  Its  tip  will  have 
at  the  tip  a  vertical  displacement 

PHj 3  Mti2 

Yi  =  -  +  -  (!) 

3EI  2EI 
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and  an  angular  displacement 


Md!  Pi!,!2 

|  0  =  -  +  -  (2) 

1  El  2EI 

I 

j  where  I  is  the  second  moment  of  area  of  the  cross-section  of  the  elastic 

j  portion  of  the  beam. 


b  i  d ! 3 

I  =  —  in4  .  (3) 

12 


Summing  the  vertical  forces  acting  on  the  rigid  portion, 

ff  “  p  '  fd  =  (4) 
and  summing  moments 

pU2-£i)-fdU3-*2)+ffU4-?,2Mi  =  J0  (5) 

where  J  is  the  second  moment  of  mass  of  the  rigid  portion  about  a  transverse 
axis  through  its  centroid.  The  mass  of  the  rigid  portion  is  m,  in  units  of 
lbf  s2in-1. 


m  =  W(c,4-a1)b2d2  x 


J 


Uir*i)2+ 


1 

32.2x12 


lb. 


1  b  ^  s2in-1 
s2  in 


(6) 

(7) 


For  small  vibrations  y2  ~y i  +  {z2-1\)b 

y  2  =  yi+  ( 1 1 )  8 
y2  =  yi+  (a2-£i  )© 
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(4)  becomes  P  =  Fp-FD-m[y,  +  (ii?-e,  )ij] 

(5)  becomes  M  =  { Fp-FD-m[y x +( «,2 - u. ,  )o] }  U2-£,) 

-F  |-J  (  Jl  3  -  £2  )'*’^r  p(?.4-^2  )_^0 

M  =  FF(a2-aI+jt„-e2  )+FD(«.l-*.2-i3+a2 ) 

-m[yi  +(^2-«  i  )o](«.2-<Li  )-J0 
m  =  FF(a4-5,1  )-FD<e3-e1  )-m[y1 +(a2-g.1)o](«.2-fc1)-Je 


Substituting  for  P  and  M  in  (1) 

FF  £, 3  Fd  £,3  me, 3 
3EI  3EI  3EI 


[y  i +( 2.2 -«-i  )e] 


Fp(£4-£a )£!2  Fd(£3-£j )£,2  m£,2 

+ - ly 

2EI  2EI  2EI 

ti2  .. 

+( £2 -2.1  )e] (£2 -2.! ) - Je 

2EI 

5  ( 2. 4  £ i  )  2 , 2  1 


-  + 

3EI 

2EI 

*i3 

_  + 

( £3-£i  )l\ 2 

3EI 

2EI 

-me!3 

me,2 

(e 

3EI  2EI 


-mj,  i  mo i 2  o.2 

+  0  -  ( £2 -2.1 ) - (£2-£,  )2 - J 

3EI  2EI  2EI 
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=  ff 


-  Fr 


+  y. 


+  0 


ill  (P 

H- 

2EI 

3 

- 

- 

il  (e 

Pi  \ 
3-  — ) 

2EI 

3 

-mPi 2 

Pi 

(P2-  — 

2EI 

3 

f-i,! 

r  / 

EI 


^•2 2  \  J 

-  —  p  i  p?  +  —  I  +  — 


6  3 


Pi 2  U 

Define  G  =  -  (8-**-  —  ) 

2EI  3 

Pi2  Hi 

H  =  -  (p3 - ) 

2EI  3 


A  = 


B  = 


-m^i  ^i 

- (P-2-  ) 

2EI  3 


EI 


f  Pi2 

2  iz  } 

,  J " 

m 

! _ _ 

-  -PjP2+  - 

1+  - 

-- 

1  6 

3  2  j 

2  . 

fd 

+  Ay  i 

+  B0  • 

Substituting  for  P  and  M  in  (2) 
F 

ci 

J  P  i 


Pi  Pi  mPi  . 

0=  fc(p4-Pi)  —  Fn(p3-Pi) - [y  +U2-Pi)e](*2“*i 

EI  u  EI  EI  1 


EI 


Fc  Pi2  Fn  Pi2  mPi 2  .. 

a  +  -E - -2 - [yi  + 


2EI 


2EI  2EI 
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(16) 

(17) 

) 

(P2-P1 )0] 
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Then  e  =  K  Fp  -  L  FQ  +  Cy;  +  De  (22) 

(17)  and  (22)  are  sufficient  to  describe  the  motion  of  the  beam. 

3.  CSMP  PROGRAM 

CSMP,  the  Continuous  System  Modeling  Program,  is  an  IBM  package 
which  performs  a  digital  simulation  of  an  analog  computer  (Reference  1). 

The  user  writes  differential  equations  to  describe  the  process  being 
simulated,  and  inputs  initial  conditions,  constants  and  other  relevant 
information  to  the  program  in  a  specified  format.  CSMP  can  list  five 
variables  as  functions  of  time,  and  can  plot  one  variable  as  a  function 
of  any  other  variable.  Run  time  for  the  CSMP  "Vibrating  Cantilever" 
program  is  typically  five  minutes,  after  which  the  user  can  change  initial 
conditions  or  constants  through  the  1130  keyboard.  The  program  demonstrates 
the  effects  of  changes  in  dimensions,  material  density  and  stiffness, 
viscous  damping  factor,  and  added  (or  subtracted)  mass  upon  the  natural 
vibration  frequencies  of  the  elastic-rigid  cantilever  beam.  It  also 
demonstrates  the  effect  of  those  changes  upon  the  transient  response  of 
the  beam  to  any  tip-applied  vertical  forcing  function. 

3.1  Equations  of  Motion 

The  equations  of  motion,  from  section  2,  are: 

yj  =  G  Fp  -  H  FD  +  Ayj  +  Be  (17) 

e  =  K  Fp  -  L  Fd  +  Cyj  +  De  (22) 

In  order  to  facilitate  construction  of  the  CSMP  block  diagram, 
the  equations  are  manipulated  as  follows: 


From  (17) 

yx  -  GFp  +  H  Fp  -  Be 
^  A 

(23) 

From  (22) 

e  -  K  Fp  +  L  FD  -  De 

yi  = 

(24) 

UNCLASSIFIED 


UNCLASSIFIED 


9. 


Therefore  Cyt  -  CGFp  +  CHF^  -  C80 
=  A0  -  KAFp  +  LAFd  -  AD0 
0(CB-AD) 

0 

Rewriting  (24), 

yi  =  C6  "  C  ®  ‘  C  FF  +  C  FD  ' 

(25)  and  (26)  are  the  equations  of  motion  presented  by  the  portion  of  the 
CSMP  block  diagram  shown  in  Figure  5. 


=  -Ae  +  (KA-CG)Fp  +  Cyj  +  (CH-LA)FQ 

,  +  _c_  v  +  |ka-cg\  p  Jch-la] 

CB-AD  0  CB-AD  yi  CB-AD  rF  \CB-AD J  *  (25) 


This  portion  of  the  CSMP  block  diagram  is  shown  in  Figure  6. 

N  is  the  damping  factor  with  units  of  lbf  s  in"1.  It  is  set  equal 
to  zero  for  undamped  vibration,  or  can  take  on  any  positive  value  to  simulate 
a  desired  degree  of  viscous  damping.  The  damper  shown  in  Figure  7  was 
assumed  to  exist  in  order  to  estimate  a  typical  magnitude  of  the  damping 
factor. 


It  was  assumed  that  there  was  a  linear  velocity  distribution  in 
the  oil  between  the  plate,  moving  with  velocity  y,  and  the  stationary 
wall.  Thus  the  velocity  gradient  was  ^  s"1.  According  to  Newton's 


law  of  viscosity 
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UNCLASSIFIED  10. 

where  t  =  shear  stress,  lbf  in"? 

u  =  viscosity  of  the  oil,  lbp  s  in-2 
=  .00142  lbf  s  in"2. 

The  total  area  of  the  plate  acted  upon  by  the  shear  stress  was 
Arlaje  =  2  x  .25  x  .25  =  .125  in2. 

The  damping  force  was 

fd  =  T  aplate 

y3 

=  .00142  x  -Q5g  x  .125 
=  .00306  x  y3  . 

FD 

From  (27) ,  N  =  — 

y3 

=  .00306  lbf  s  in"1  . 

3. 3  Forcing  Function 

The  force  applied  to  the  tip  of  the  beam,  Fp,  was  modelled  as  a 
linear  rise  to  the  peak  force,  Fp^*  at  time  t  =  .IT,  followed  by  an 
exponential  decay.  Figure  8  illustrates  the  nondimensionalized  forcing 
function  Fp/FFMAX  versus  t/T. 


From  t/T  =  0  to  t/T=. 1 , 

!l_  =  10 1 

F  u  T 

rFMAX  1 


From  t/T  =  .1  to  t/T  1.0, 


*  ru  n 
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I 

! 


t> 


To  evaluate  the  constant  k  it  was  assumed  that 

Fp/FfMAx  =  ‘2  when  t//T  = 

.  2  =e-k(-6-^ 

=  e-5k 

k  = 

-.5 

=  3.219  . 


Therefore 


.  e-3‘219  f  "  -1 


FMAX 


,  .1  <  j  <  l.o  . 


Table  II  lists  values  of  Fp/FFMAX  at  regular  intervals  of  t/T. 


TABLE  II 


-  vs 

hFMAX 


FFMAX 


0. 

0. 

.1 

1.0 

.2 

.7246 

.3 

.5252 

.4 

.3806 

.5 

.2758 

.6 

.2000 

.7 

.1448 

.8 

.1050 

.9 

.0760 

1.0 

.0552 
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The  forcing  function  part  of  the  CSMP  block  diagram  is  shown 
in  Figure  9.  In  the  CSMP  program,  the  values  of  Table  2  are  read  in  as 
specifications  for  a  Function  Generator  Block,  Block  17.  The  times  of 
the  beginning  and  end  of  the  force  pulse,  i.e.,  0  s  and  T  s,  must  also 
be  specified  as  Parameters  2  and  1,  respectively,  of  the  Function  Gene¬ 
rator  Block.  Input  to  the  Function  Generator  is  t  from  the  time  base. 
Block  76.  Output  is  Fp/FpM^,.  When  t/T  is  not  a  multiple  of  .1,  the 

Function  Generator  linearly  interpolates  a  value  for  Fp/Fp^.  For  t/T 

greater  than  1.0,  Fp/Fp^x  is  constant  at  .0552.  Fp^  is  specified  as 

Parameter  1  of  a  Gain  Block,  Block  18,  which  receives  Fp/Fp^  as  input, 
and  which  outputs  Fp. 

3.4  Strain 


It  was  desired  that  the  CSMP  program  should  continuously  cal¬ 
culate  the  maximum  strain  at  the  base  of  the  beam,  because  strain  gauges 
were  placed  there  on  the  actual  beam.  Linear  elastic  deformation  was 
assumed,  and  the  dead  weight  of  the  beam  was  neglected. 


_  o 
e  -  E 

where  e  =  strain,  in  in-1 
o  =  stress,  lbf  in-2 

E  =  Young's  modulus  of  beam  material,  lb^.  in-2. 
Applying  the  basic  flexure  formula, 


(28) 


E  E  x 


M0  “i 

~TT 


(29) 


where  M0  =  bending  moment  at  base  of  beam,  lbp  in 
dj  =  thickness  of  beam  at  its  base,  in 
I  =  second  moment  of  area  of  beam  cross-section  at 
its  base,  in*4.  See  equation  (3). 

From  the  free  body  diagram.  Figure  4,  it  can  be  shown  that 

M0  =  Pt,  +  M  .  (30) 


- - .i ■  .  TW— 
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From  equation  (5), 

M  =  PU2-*i)  -  Fq( a3-2.2)  +  FfU4-22)  -  Jo  .  (31) 

Recalling  equation  (11), 

P  =  Fp  -  FD  -  m[yj  +  (fc2-n1)e]  .  (11) 

Equating  the  left-hand  side  of  (29)  to  the  right-hand  side  of 
(30)  and  substituting  for  M0,  P  and  M, 


e  - 


j-|"m|ff  -  fd  “  m[yi  +  U2“*i)0]  j  + 

P(fi-2-fi-l )  -  FD(^3-fc2)  -  PpK'^)  -  06  J  • 


(32) 


Equation  (32)  is  represented  by  the  portion  of  the  CSMP  block 
diagram  shown  in  Figure  10. 

Figure  11  shows  the  complete  CSMP  block  diagram.  Appendix  A  is 
a  listing  of  the  CSMP  program  for  a  typical  run. 


4.  FORTRAN  PROGRAM 

The  main  purpose  of  the  FORTRAN  program  is  to  punch  Initial 
Condition/Parameter  cards  for  the  CSMP  program.  It  also  prints  a  listing 
of  the  parameters.  Only  the  parameters  dependent  on  beam  dimensions  are 
calculated  by  this  program. 

Figure  12  shows  a  flow  chart  of  the  FORTRAN  program.  It  begins 
by  reading  a  beam  dimension  data  card.  Parameters  needed  for  the  CSMP 
simulation  of  this  particular  beam  are  calculated  and  stored  in  an  output 
matrix.  The  program  reads  another  card,  calculates  and  stores  the  next 
set  of  parameters,  and  so  on  until  a  blank  card  is  read.  Then  the  para¬ 
meter  listing  is  printed  and  the  Initial  Condition/Parameter  cards  for  the 
CSMP  program  are  punched. 

Table  III  defines  the  variables,  and  Appendix  B  is  a  listing  of 
the  program  with  a  sample  of  output. 
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The  variable  ADMAS  represents  mass  added  to  the  beam  at  the 
centre  of  gravity  of  the  rigid  portion.  A  nut  and  bolt  or  a  lump  of 
solder  might  be  attached  to  the  beam  to  lower  the  natural  frequencies  of 
vibration.  The  programs  are  only  valid  for  mass  added  at  the  C  of  6 
of  the  rigid  portion  because  they  neglect  any  change  in  the  second  moment 
of  mass,  J,  due  to  the  added  mass. 

TABLE  III 

FORTRAN  VARIABLES 


Variable 

Definition 

Mode 

Class 

U(K,N) 

beam  dimension  input  matrix 

Real 

2 

V(L,N) 

parameter  output  matrix 

II 

2 

EMM(N) 

beam  mass,  g 

II 

1 

N 

index  to  distinguish  data  sets 

Integer 

0 

ELI 

*i>  in 

Real 

0 

EL2 

ii2,  in 

«l 

0 

EL3 

«.3,  in 

II 

0 

EL4 

8-4,  in 

II 

0 

B1 

bi,  in 

II 

0 

B2 

b2,  in 

II 

0 

Dl 

di,  in 

II 

0 

D2 

d2,  in 

II 

0 

W 

w,  lbf  in'3 

II 

0 

E 

E,  lbj.  in-2 

ll 

0 

ADMAS 

mass  added  to  C  of  G  of  rigid 

portion  of  beam,  g 

II 

0 

EM 

m,  lbf  s2in_1.  See  equation  (6) 

II 

0 

REALJ 

J,  lb^  s2in.  See  equation  (7) 

II 

0 

REAL  I 

I,  in4.  See  equation  (3) 

II 

0 

A 

A,  s2.  See  equation  (15) 

II 

0 

B 

B,  in  s2.  See  equation  (16) 

II 

0 
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TABLE  III  (continued) 


Variable 

Defini tion 

Mode 

Class 

C 

C,  in-1  s2.  See  equation  (20) 

li 

0 

D 

D,  s2.  See  equation  (21) 

li 

0 

G 

G,  in  lb^-1.  See  equation  (13) 

II 

0 

H 

H,  in  lb^-1.  See  equation  (14) 

li 

0 

REALK 

K,  lb^-1.  See  equation  (18) 

II 

0 

REALL 

L,  lb^.-1.  See  equation  (19) 

li 

0 

P1B3 

parameter  1,  block  3,  s2 

Real 

0 

P1B4 

parameter  1,  block  4,  in 

li 

0 

P2B5 

parameter  2,  block  5,  s2 

li 

0 

P1B6 

parameter  1,  block  6,  in  s-2 

li 

0 

P2B6 

parameter  2,  block  6,  Ib^*1  in  s~ 

2  ■' 

0 

P3B6 

parameter  3,  block  6,  Ib^-1  in  s" 

2  ■' 

0 

P1B7 

parameter  1,  block  7,  lbf  s2  in-1 

ii 

0 

P1B10 

parameter  1,  block  10,  in 

it 

0 

P2B10 

parameter  2,  block  10,  in 

» 

0 

P3B10 

parameter  3,  block  10,  in 

ii 

0 

P1B11 

parameter  1,  block  11,  in 

n 

0 

P1B12 

parameter  1,  block  12,  in"1  lbf-1 

ii 

0 

P2B13 

parameter  2,  block  13,  in 

ii 

0 

P1B15 

parameter  1,  block  15,  in-1  s"2 

ii 

0 

P2B15 

parameter  2,  block  15,  lbf_1  s"2 

ii 

0 

P3B15 

parameter  3,  block  15,  lb^"1  s-2 

n 

0 

P2B16 

parameter  2,  block  16,  lb^  s2  in 

ii 

0 

I  total  number  of  data  sets  Integer  0 
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5.  USE  OF  PROGRAMS 
5.1  CSMP 

Reference  1  gives  a  detailed  description  of  how  to  run  a  CSMP 
program.  Only  a  brief  outline  is  given  here. 

1.  Prepare  a  card  deck  which  includes  appropriate  control 

cards:  Configuration  Specification  cards.  Initial 
Condition/Parameter  Specification  cards.  Function 
Generator  Specification  cards  and  about  a  dozen  blank 
cards. 

2.  Ensure  that  the  correct  disc  is  installed. 

3.  Load  the  card  deck  and  start  the  program. 

4.  The  console  printer  will  list  the  Configuration 

Specifications,  I .C. /Parameter  Specifications  and 
Function  Generator  Specifications.  It  will  list  the 
console  switches  to  be  operated,  after  a  run,  when 
choosing  program  options. 

5.  The  console  printer  will  ask  the  operator  to  enter 

values  for: 

Integration  Interval 
Total  Time 

Block  for  Y-Axis  of  Plot 

Minimum  and  Maximum  Values  on  Y-Axis 

Block  for  X-Axis  of  Plot 

Minimum  and  Maximum  Values  on  X-Axis. 

6.  Press  the  "Start"  button  on  the  console  and  the 

plotter  will  draw  the  plot  frame. 

7.  The  console  printer  will  ask  the  operator  to  enter: 

Print  Interval 

Block  Numbers  for  five  variables  to  be  listed 
on  the  line  printer. 

As  soon  as  the  last  Block  Number  is  entered,  the  run 
will  automatically  begin. 
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8.  The  run  will  continue  until  the  "Total  Time"  is  reached, 

or  may  be  stopped  sooner  by  operating  Switch  0  on  the 
console. 

9.  If  another  run  is  required,  operate  console  switches  to 

choose  program  options.  For  example,  a  new  plot  frame 
will  be  drawn  if  Switch  7  is  operated,  and  Initial 
Conditions/Parameters  can  be  changed  if  Switch  2  is 
operated.  Then  press  "Start". 

10.  If  Configuration,  I ,C. /Parameter  or  Function  Generator 

Specifications  are  to  be  changed,  the  program  will 
attempt  to  read  new  specifications  from  the  card  deck. 
After  a  blank  card  is  read,  the  program  will  ask  for 
input  from  the  console. 

Other  specifications  are  changed  through  the  console  only. 

11.  Press  "Interrupt  Request"  to  stop  the  program  after  the 

last  run. 

Table  IV  shows  how  several  runs  could  be  set  up  to  predict 
the  behaviour  of  a  given  beam  under  various  conditions. 

Selection  of  initial  y}  and  e  is  arbitrary.  However,  a  conmon 
method  of  observing  the  free  vibration  of  an  actual  beam  would  be  to  deflect 
the  tip  a  certain  amount,  y4,  then  suddenly  release  it.  Appendix  C  explains 
how  to  calculate  initial  yY  and  o  when  the  initial  tip  deflection,  y4,  is 
specified. 

Choice  of  Integration  Interval  involves  a  compromise  between  run 
time  and  the  accuracy  of  the  numerical  integration  process  performed  by 
the  CSMP  program.  The  best  interval  is  the  longest  one  for  which  the 
integration  process  is  stable.  If  the  interval  is  too  long,  the  amplitude 
of  the  predicted  vibrations  increases  with  time  even  when  no  external 
force  is  applied,  an  obvious  error.  Cases  with  no  damping  require  the 
shortest  Integration  Intervals  because  the  frequencies  of  vibration  are 
highest.  A  good  first  guess  for  the  Integration  Interval  is  one-tenth 
the  period  of  the  highest  frequency  vibration  expected. 
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5.2  FORTRAN 

1.  Prepare  a  card  deck  as  listed  in  Appendix  B,  with  beam 

data  cards  in  the  format  illustrated  by  Table  V.  This 
program  accepts  up  to  five  beam  data  cards.  Change  the 
DIMENSION  statement  if  more  cards  are  to  be  processed 
in  one  run.  Following  the  n  beam  data  cards  should  be 
at  least  lln  blanks  on  which  will  be  punched  the  CSMP 
I .C. /Parameter  Specifications. 

2.  Load  card  deck  and  start  the  program. 

3.  When  the  run  is  finished,  transfer  the  newly-punched 

I .C. /Parameter  cards  to  the  CSMP  deck. 

6.  RESULTS 

The  CSMP  program  was  used  to  predict  frequencies  of  undamped 
free  vibration  for  the  Mark  1  and  Mark  2  beams.  Table  VI  lists  some  of 
the  parameters,  while  Figures  13  and  14  are  strain  vs  ';ime  plots.  As 
expected  from  a  two-degree-of-freedom  system,  there  are  two  waveforms 
superimposed.  Good  agreement  between  predicted  and  experimental  fre¬ 
quencies  is  shown  in  Table  VI.  The  higher  frequencies  were  not  observed 
experimentally  due  to  the  resolution  of  the  oscilloscope  on  which  strain 
vs  time  was  displayed. 

A  third  CSMP  run  predicted  frequencies  of  undamped  free  vibration 
for  a  Mark  1  beam  with  a  .500  g  mass  added  at  the  centre  of  gravity  of  the 
rigid  portion.  The  results  for  this  run  are  shown  in  Table  VI  and  Figure 
15. 


7.  CONCLUSIONS 

The  CSMP  "Vibrating  Cantilever"  program  made  good  predictions 
of  the  undamped  free  vibration  frequencies  of  the  Mark  1  and  Mark  2  beams. 
It  can  therefore  be  used  with  confidence  to  make  similar  predictions  about 
other  elastic-rigid  cantilevers. 

Further  verification  of  the  program  might  involve  the  comparison 
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of  predicted  to  experimental  strain  amplitudes  during  undamped  free  vibra¬ 
tion.  The  most  rigorous  test  would  be  a  comparison  of  predicted  to  experi¬ 
mental  damped  transient  response  to  the  forcing  function.  Accurate  esti¬ 
mates  of  damping  factor,  peak  force  and  duration  of  the  force  pulse  would 
have  to  be  supplied  to  the  program. 

Leaving  aside  the  necessity  for  verification,  the  program  can  be 
used  to  predict  the  effect  on  unwanted  oscillations  of  changes  in: 

1)  beam  dimensions 

2)  mass  added  or  subtracted  at  the  centre  of  gravity 
of  the  rigid  (flanged)  section 

3)  viscous  damping  factor 

4)  rate  of  loading 
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TABLE  VI 

UNDAMPED  FREE  VIBRATION 


Beam 

Dimensions 

Run  1 

Mark  1  Beam 

Run  2 

Mark  2  Beam 

Run  3 

Mark  1  + 

.500  g  Mass 

Beam  Dimensions 

Hi  in 

.750 

.600 

.750 

*2  " 

1.875 

1.425 

1.875 

*3  " 

2.500 

1.750 

2.500 

*4  " 

3.000 

2.250 

3.000 

bi  ” 

.250 

0.250 

.250 

b2  " 

.375 

.500 

.375 

di  11 

.010 

.020 

.010 

d2  " 

.010 

.020 

.010 

Values  assumed 

for  jw  lbf  in'3 

.295 

.295 

.295 

phosphor  bronze 

\E  1 b^  in-2 

16100000. 

16100000. 

16100000. 

ADMAS  g 

0. 

0. 

.500 

Integration  Interval  s 

.00005 

. 00001 

. 00005 

Total  Time 

s 

.050 

.050 

.050 

Initial  yk 

in 

.025 

.025 

.025 

y\ 

in 

. 00343 

.00525 

. 00343 

"  9 

rad 

.00938 

.01197 

.00938 

N 

lbf  s  in-1 

0. 

0. 

0. 

ffmax 

lbf 

0. 

0. 

0. 

Predicted  Frequencies  of 

Vibration 

Hz 

25 

74 

22 

550 

1500 

500 

Experimental  Frequencies  of 

Vibration 

Hz 

28 

60  to  70 

“  — 
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Figure  2:  Mark  1  Beam  and  Simplified  Model 
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,  =  Vertical  displacement  of  elastic-rigid  junction  from  rest  position,  in 

-  Vertical  displacement  of  C  of  6  of  rigid  position  from  rest  position,  in 

y,  =  Vertical  displacement  of  damping  force  point  of  application  from 
rest  position,  in 

y i,  =  Vertical  displacement  of  tip 

e  =  Angular  displacement  of  rigid  position  of  beam  from  rest  position,  rad 


Figure  3:  Coordinates 
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NOTES:  1.  Arrowheads  represent  connections  with  other  parts  of  the 
block  diagram. 

2.  Blocks  20,  21,  22,  23  and  24  are  gain  factors  which  allow 
a  sufficient  number  of  significant  digits  when  reading  in 
the  parameters  for  blocks  3,  5  and  15. 


Figure  5:  CSMP  Block  Diagram  -  Equations  of  Motion 
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3:  CSMP  Block  Diagram  -  Viscous  Damping 
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Figure  9:  CSMP  Block  Diagram  -  Forcing  Function 
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1.  Blocks  25  and  26  are  gain  factors  which  allow  a  sufficient 
number  of  significant  digits  when  reading  in  the  parameters 
for  blocks  16  and  7. 


Figure  10:  CSMP  Block  Diagram  -  Strain 
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FIG.  13 :  STRAIN  VS  TltfE-  RUN  I 
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FIG.  14:  STRAIN  VS  TIME  -  RUN  2 
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FIG.  I5:  STRAIN  VS  TIME  -  RUN  3 
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FUNCTION  GENERATOR  SPECIFICATIONS 

0.0000  1.0000  0.7246  0.5252 
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APPENDIX  B  -  Fortran  Program  Listing 
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Summing  moments  on  the  rigid  portion  of  the  beam, 

Mi  =  .  {2 

The  deflection  and  rotation  at  the  end  of  the  elastic  portion 

Fi  M3  Fi(«4-^1)^12 


=  Fi  ^  !  3  (8-4-fi.x  )lti: 
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Assuming  that  the  rotation  0  is  small, 

yi  =  y4  -  U4-*i)0 
Substituting  equation  (36)  in  (34) 
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Rearranging  (37) 


i  Pil 

4_^i )  L  3 


(*4-*i)*i2 


Equating  the  right-hand  side  of  (35)  and  (38) 
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